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We study the kicked rotator in the classically fully chaotic regime using Izrailev's A^'-dimensional 
model for various A'^ < 4000, which in the limit N ^ oo tends to the quantized kicked rotator. We 
do not treat only the case K = 5 as studied previously, but many different values of the classical 
kick parameter 5 < K < 35, and also many different values of the quantum parameter k £ [5, 60]. 
We describe the features of dynamical localization of chaotic eigenstates as a paradigm for other 
both time-periodic and time-independent (autonomous) fully chaotic or/and mixed type Hamilton 
systems. We generalize the scaling variable A = loo/N to the case of anomalous diffusion in the 
classical phase space, by deriving the localization length loo for the case of generalized classical 
diffusion. We greatly improve the accuracy and statistical significance of the numerical calculations, 
giving rise to the following conclusions: (CI) The level spacing distribution of the eigenphases (or 
quasienergies) is very well described by the Brody distribution, systematically better than by other 
proposed models, for various Brody exponents /3br. (C2) We study the eigenfunctions of the Floquet 
operator and characterize their localization properties using the information entropy measure, which 
after normalization is given by /3ioc in the interval [0, 1]. The level repulsion parameters /3br and /3ioc 
are almost linearly related, close to the identity line. (C3) We show the existence of a scaling law 
between Pioc and the relative localization length A, now including the regimes of anomalous diffusion. 
The above findings are important also for chaotic eigenstates in time-independent systems (Batistic 
and Robnik 2010,2013), where the Brody distribution is confirmed to a very high degree of precision 
for dynamically localized chaotic eigenstates even in the mixed-type systems (after separation of 
regular and chaotic eigenstates). 
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I. INTRODUCTION 

One of the main cornerstones in the development of 
quantum chaos [H-Q is the finding that in classically fully 
chaotic, ergodic, autonomous Hamilton systems vi^ith the 
purely discrete spectrum the fluctuations of the energy 
spectrum around its mean behaviour obey the statistical 
laws described by the Gaussian Random Matrix The- 
ory (RMT) [3, [1] , provided that we are in the sufficiently 
deep semiclassical limit. The latter condition means that 
all relevant classical transport times are smaller than 
the so-called Heisenberg time, or break time, given by 
tn — 2Tih/ IS.E, where h = 27rfi, is the Planck constant and 
Ai? is the mean energy level spacing, such that the mean 
energy level density is p{E) = 1/AE. This statement is 
known as the Bohigas - Giannoni - Schmit (BGS) con- 
jecture and goes back to their pioneering paper in 1984 
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d, although some preliminary ideas were published in 
31 . Since AE esc , where / is the number of degrees of 
freedom (= the dimension of the configuration space), we 
see that for sufficiently small h the stated condition will 
always be satisfied. Alternatively, fixing the h, we can go 
to high energies such that the classical transport times 
become smaller than tn- The role of the antiunitary 
symmetries that classify the statistics in terms of GOE, 
GUE or GSE (ensembles of RMT) has been elucidated in 
@] , see also [9] , and [J-Ill . The theoretical foundation for 



the BGS conjecture has been initiated first by Berry 1C|, 
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and later further developed by Richter and Sieber 
arriving finally in the almost-final proof proposed by the 
group of F. Haake ^^t^- 

On the other hand, if the system is classically inte- 
grable, Poisson statistics applies, as is well known and 
goes back to the work by Berry and Tabor in 1977 (see 
[l|-|3| and the references therein, and for the recent ad- 
vances [T6j). 

In the mixed type regime, where classical regular 
regions coexist in the classical phase space with the 
chaotic regions, being a typical KAM-scenario which is 



2 



the generic situation, the so-cahed Principle of Uniform 
Semiclassical Condensation (of the Wigner functions of 
the eigenstates; PUSC) apphes, based on the ideas by 
Berry [ll], and further extended by Robnik Conse- 
quently the Berry- Robnik statistics JJi, J^] is observed - 
see also [3| - again under the same semiclassical condition 
stated above requiring that tn is larger than all classical 
transport times. 

The relevant papers dealing with the mixed type 
regime after the work fl^ are [20l - [27j and the most recent 
advance was published in [2^, while [l^ is the relevant 
work in progress. If the couplings between the regular 
eigenstates and chaotic eigenstates become important, 
due to the dynamical tunneling, we can use the ensembles 
of random matrices that capture these effects • the 
tunneling strengths typically decrease exponentially with 
the inverse effective Planck constant, they rapidly disap- 
pear with increasing energy, or by decreasing the value 
of the Planck constant. 

Here it must be emphasized that the analogies between 
the time-periodic systems (the kicked rotator) and time- 
independent systems (like mixed-type billiards) that we 
are drawing and studying in this paper refer to the 
chaotic eigenstates only, which means that we have to 
conceptually separate the regular and the chaotic eigen- 
states in each system. If the semiclassical condition is 
satisfied, then for the subspectrum of the chaotic eigen- 
states we find extendedness and GOE statistics. This 
should be compared with the extended states in finite 
dimensional kicked rotator model for K > 7, where the 
corresponding classical dynamics is fully chaotic. 

However, if the semiclassical condition is not satisfied, 
such that tn is no longer larger than the relevant classical 
transport time, like e.g. the diffusion time in fully chaotic 
but slowly ergodic systems, we find the so-called dynam- 
ical localization (or Chirikov localization) first ob- 
served in time-dependent systems (see e.g. [ij), which are 
the main topics of the present work and will be discussed 
below in detail, but later on analyzed quite systematically 
in autonomous (time-independent) systems by many au- 
thors. For an excellent review see the paper by Prosen 
3l| and the references therein. In such a situation it 
turns out that the Wigner functions of the chaotic eigen- 
states no longer uniformly occupy the entire classically 
accessible chaotic region in the classical phase space, but 
are localized on a proper subset of it. In contradistinc- 
tion to the tunneling effects, these dynamical localization 
effects can survive to very high lying eigenstates. Indeed, 
this has been analyzed with unprecedented precision and 
statistical significance by Batistic and Robnik 28] in case 
of mixed type systems, and the work is being extended in 
the analysis of separated regular and chaotic eigenstates 
[29, 32]. The most important discovery is that the level 
spacing distribution of the dynamically localized chaotic 
eigenstates is very well described by the Brody distri- 
bution, introduced in [s^, see also [s^l, with the Brody 
parameter values /3br within the interval [0,1], where 
/3br = yields the Poisson distribution in case of the 



strongest localization, and /3br = 1 gives the Wigner sur- 
mise (2D GOE, as an excellent approximation of the infi- 
nite dimensional GOE) . To our great surprise the Brody 
distribution fits the empirical data much better than the 
distribution function proposed by F. Izrailev (see [35l.l36j 
and the references therein) characterized by the param- 
eter /3iz. This is still true also for the improved Izrailev 
distribution published in [37] and recently used in 38| . It 



is well known that Brody distribution so far has no the- 
oretical foundation, but our empirical results show that 
we have to consider it seriously thereby being motivated 
for seeking its physical foundation. 

In the present study of the kicked rotator, besides the 
above mentioned results on the relevance of the Brody 
distribution, we go beyond Izrailev's results in that we 
study not only the case of the classical kick parameter 
K = 5, but for many other K g [5,35], and many dif- 
ferent values of the quantum parameter fc, and consider 
the relevance of the classical diffusion in greater depth, 
allowing also for the anomalous diffusion. In so doing we 
greatly generalize and improve the evidence for the linear 
relationship between the information entropy localization 
measure /3ioc and /3br, and also for the scaling relation- 
ship between /3ioc and the scaling variable A, which is the 
theoretical localization length divided by the dimension 
of the system. 

Our work corroborates the view (see [s^ and the refer- 
ences therein, especially (39l - l4l| ) that time-independent 
and time-periodic chaotic systems have many properties 
in common when it comes to the statistical properties 
of discrete energy spectra and the discrete quasienergy 
spectra (or eigenphases), respectively. We think that this 
view can be extended also to quantifying the degree of 
localization in such systems. Very recent results (Batistic 
and Robnik 2013) confirm this expectation, and will be 
published separately. 

The paper is organized as follows: In section |ll] we 
introduce the model system, in section IIIII we describe 
the aspects of generalized diffusion in the classical sys- 
tem (standard map) including the accelerator modes and 
the anomalous diffusion and relate it to the quantum lo- 
calization properties, deriving the new formula for the 
localization length. In section IIVI we define the finite di- 
mensional quantum model system, introduced by Izrailev 
and study not only the cases of quantum resonance, but 
also the generic cases. In section|V]we define the informa- 
tion entropy localization measure Aoc, in section IVll we 
study the statistical properties of spectra (eigenphases), 
in section IVIII we analyze the relationship between the 
localization parameter /3ioc and the spectral level repul- 
sion parameters /3br and /3iz, and also study the scaling 
relationship between /3ioc and the scaling parameter A. 
In section IVIIII we draw the final conclusions and discuss 
the results. In Appendix [X] we define and explain the 
[/-function of the level spacings and in Appendix |B] we 
show some additional results on energy level statistics of 
the quantum kicked rotator. 
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II. INTRODUCING THE MODEL 

One of the main models of time-dependent systems is 
the kicked rotator introduced by Casati, Chirikov, Ford 
and Izrailev in 1979 We introduce it here in detail 
for the purpose of defining and fixing the variables and 
the notation. The Hamiltonian function is 



H=^ + VoSTit) cos9. 



(1) 



It is one of the most important paradigms of classical 
conservative (Hamiltonian) systems in nonlinear dynam- 
ics. Here p is the (angular) momentum, / the moment 
of inertia, Vq is the strength of the periodic kicking, 9 is 
the (canonically conjugate, rotation) angle, and Sxit) is 
the periodic Dirac delta function with period T. Since 
between the kicks the rotation is free, the Hamilton equa- 
tions of motion can be immediately integrated, and thus 
the dynamics can be reduced to the standard mapping, 
or so-called Chirikov- Taylor mapping, given by 



Pn+l = Pn + Vo sin6'„+i, 9n+l = 



T 

jPn, 



(2) 



and introduced in |43l - l45l | . Here the quantities {9n,Pn) 
refer to their values just immediately after the n-th kick. 
Obviously, by introducing new dimensionless momentum 
Pn = PnT/I, we get 



Pn+l = Pn +Ksin9n+l, 9n+l = 9„ + Pn, 



(3) 



where the system is now governed by a single classical 
dimensionless control parameter K = VqT/I, and the 
mapping is area preserving. 

The quantum kicked rotator (QKR) is the quantized 
version of Eq. ([1]) , namely 



^2 g2 

^ = -^5^ + ^o5t(0 cosl 



(4) 



The physics of the QKR is extremely rich and it is a 
paradigm of quantum chaos in Floquet {— time-periodic) 
systems 'sg']. It is also relevant for the autonomous 
Hamilton systems as indicated above. For such a Floquet 
system the Floquet operator F acting on the wavefunc- 
tions (probability amplitudes) ip{9), 9 G [0,27r), upon 
each period (of length T) can be written as (see e.g. [l|. 
Chapter 4) 



F 



exp 



iVo 



■ COS! 



exp 



ihT 



(5) 



where now we have obviously two dimensionless quantum 
control parameters 

TiT 



k = 



(6) 



which satisfy the relationship K = kr = VqT/I, K being 
the classical dimensionless control parameter of Eq. ([3|) . 
By using the angular momentum eigenfunctions 



\n) = an{9) 



exp(i n 9), 



(7) 



where n is any integer, we find the matrix elements of F , 
namely 



Fmn^ {m\F\n) = exp 



-n^ 1 z"-™J„-™(fc), (8) 



where Jv{k) is the v-th order Bessel function. For a 
wavefunction '0(^) shall denote its angular momen- 
tum component (Fourier component) by 



2tv 



<{0)md9 



27r 7o 



2Tr 



i){9) exp{-in9) d9. 



(9) 



The QKR has very complex dynamics and spectral prop- 
erties. As the phase space is infinite (cylinder), p g 
(— oo,+oo),6' e [0,27r), the spectrum of the eigenphases 
of F, denoted by (/)„, or the associated quasienergies 
hwn — fitpn/T, introduced by Zeldovich [i^, can be con- 
tinuous, or discrete. It is quite well understood that for 
the resonant values of r 



47rr 



(10) 



q and r being positive integers without common fac- 
tor, the spectrum is continuous, as rigorously proven by 
Izrailev and Shepelyansky |47h50I|. and the dynamics is 
(asymptotically) ballistic, meaning that starting from an 
arbitrary initial state the mean value of the momentum 
(p) increases linearly in time, and the energy of the sys- 
tem E = (p^)/(2/) grows quadratically without limits. 
For the special case q = r = 1 this can be shown elemen- 
tary. Such behaviour is a purely quantum effect, called 
the quantum resonance. Also, the regime of quadratic 
energy growth manifests itself only after very large time, 
which grows very fast with the value of the integer q from 
Eq. pn]) . such that for larger q this regime practically 
cannot be observed. 

For generic values of r/(47r), being irrational number, 
the spectrum is expected to be discrete but infinite. But 
the picture is very complicated. Casati and Guarneri 
5l| have proven, that for r/(47r) sufficiently close to a 
rational number, there exists a continuous component 
in the quasienergy spectrum. So, the absence of dy- 
namical localization for such cases is expected as well. 
Without a rigorous proof, we finally believe that for all 
other ( "good" ) irrational values of r/ (Att) we indeed have 
discrete spectrum and quantum dynamical localization. 
In such case the quantum dynamics is almost periodic, 
and because of the effective finiteness of the relevant set 
of components m„ and of the basis functions involved, 
just due to the exponential localization (see below), it is 
even effectively quasiperiodic (effectively there is a finite 
number of frequencies), and any initial state returns after 
some recurrence time arbitrarily close to the initial state. 
Thus the energy cannot grow indefinitely. 
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III. LOCALIZATION AND DIFFUSION 
PROPERTIES 

In the generic (nonresonant) case we observe in the 
semiclassical regime of large fc ^ 1 and in the classi- 
cally chaotic regime K > Kcrit ~ 0.9716..., the so-called 
dynamical localization also called Chirikov localization: 
starting from an initial semiclassical wave packet of the 
width smaller than the localization length, to be precisely 
defined below, the average energy grows first according 
to the classical diffusion, but stops after a finite time, 
i.e. the localization time iioc (physical time divided by 
the period of kicking T, that is the number of kicks, and 
thus dimensionless), which is derived below. 

The asymptotic localized eigenstates are quasistation- 
ary, they just oscillate under the action of F, as the quan- 
tum recurrence time is very large. They are very special, 
as their expansion coefficients in the basis of the angu- 
lar momentum eigenstates \n) must be highly correlated. 
In fact, more can be said about these asymptotic eigen- 
states: they are exponentially localized. The (dimension- 
less) localization length in the space of the angular mo- 
mentum quantum numbers is derived below, and is equal 
(after introducing some numerical correction factor a^) 
to the dimensionless localization time t\oc [Eq. p7|) ] . We 
denote it like in reference [36j by lea- Therefore, an ex- 
ponentially localized eigenfunction centered at m in the 
angular momentum space [Eq. ([7])] has the following form 



1 



In 



■ exp 



2|m- 



(11) 



where u„ is the probability amplitude [Eq. of the 
localized wavefunction ipiO). The argument leading to 
tioc in Eq. ([T7)) originates from the observation of the 
dynamical localization by Casati et al j42j . and in par- 
ticular from [s^ , and and is well explained in [l[ , in case 
of normal diffusion = 1 , whilst for general fi we give a 
theoretical argument in this section. 

Since the spectrum is discrete we can ask questions 
about the statistical properties of the spectrum of the 
quasienergies, or eigenphases. However, since the sys- 
tem is infinite with infinitely many exponentially local- 
ized eigenstates, we have infinitely many eigenphases on 
the interval [0,27r), resulting in an infinite level density, 
and thus all level spacings are zero. Nevertheless, for any 
finite but arbitrarily large number of eigenstates N ev- 
erything is well defined. In the classically fully chaotic 
regime one would naively expect the applicability of the 
RMT, in our case the GOE statistics. However, this is 
not observed. On the contrary, the statistics is Poisso- 
nian, which is the consequence of the finite localization 
length at any k in the infinite basis of the angular mo- 
mentum. Following the heuristic argument by Izrailev, 
we can say that eigenstates can be quasienergetically very 
close to each other, an almost degenerate pair, although 
they are located in the angular momentum space very 
far from each other and practically do not overlap due 



to the exponential localization. Therefore, they do not 
"feel" each other, they do not interact, in the sense that 
changing slightly one of them does not change the other 
one, and thus contribute to the spectrum in a completely 
uncorrelated way. This results in the Poissonian statis- 
tics. 

The question arises, where do we see the analogous 
phenomena predicted by the RMT and observed in the 
quantum chaos of time-independent bound systems with 
discrete spectrum. To see these effects the system must 
have effectively finite dimension. Truncation of the in- 
finite matrix Fmn in Eq. ([5]) in tour de force is not ac- 
ceptable, even in the technical case of numerical compu- 
tations, since after truncation the Floquet operator is no 
longer unitary. 

The only way to obtain a quantum system which shall 
in this sense correspond to the classical dynamical sys- 
tem [Eqs. (P), ^ and ©] is to introduce a finite iV- 
dimensional matrix, which is symmetric unitary, and 
which in the limit N ^ oo becomes the infinite dimen- 
sional system with the Floquet operator [Eq. ([S])]. The 
semiclassical limit is fc — oo and r — > 0, such that 
K ~ kr = constant. As it is well known [36|, for the 
reasons discussed above, the system behaves very simi- 
larly for rational and irrational values of r/ (47r). Such a 
A'^-dimensional model [35| will be introduced in the next 
section IIVI 

Let us now derive the semiclassical estimate of the lo- 
calization time t\oc and the localization length /qo, both 
being dimensionless. It turns out that they are equal, as 
shown in Eq. (jl6|) . The generalized diffusion process of 
the standard map [Eq. ^] is defined by 

{{APr)^D,{K)n'^, (12) 

where n is the number of iterations (kicks) , and the ex- 
ponent /i is in the interval [0,2), and all variables P, 9 
and K are dimensionless. Here D^{K) is the general- 
ized classical diffusion constant. In case fj, ~ 1 we 
have the normal diffusion, and Di (K) is then the normal 
diffusion constant, whilst in case of anomalous diffusion 
we observe subdiffusion when < < 1 or superdiffusion 
if 1 < /i < 2. In case /i = 2 we have the ballistic trans- 
port which is associated with the presence of accelerator 
modes (see below). 

As the real physical angular momentum p and P are 
connected by P = pT/I we have for the variance of p the 
following equation 



{{Apr) = ^D,r 



(13) 



Now we argue as follows: The general wisdom (golden 
rule) in quantum chaos is that the quantum diffusion 
follows the classical diffusion up to the Heisenberg time 
(or break time, or localization time), defined as 

tH = ^, (14) 

where AE is the mean energy level spacing. In our 
case we have the quasienergies and AE = hAuj, where 
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Aw = A(/)/T, and A(j> is the mean spacing of the eigen- 
phases. This might be estimated at the first sight as 
A(j) — 271 /N, but this is an underestimate, as effectively 
we shall have due to the locaUzation only Zoo levels on 
the interval [0, 2tt). Therefore Ac/) — 2tt/Ioo and we find 



tH 



2'kT 



TI 

J- Ir: 



(15) 



Since T is the period of kicking, and tn is the real phys- 
ical continuous time, we get the result that the discrete 
time (number of iterations of Eq. ^ at which the quan- 
tum diffusion stops), the localization time iioc is indeed 
equal to the localization length in momentum space, i.e. 
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Since our derivation is not rigorous, we use the approx- 
imation symbol rather than equality, in particular as 
the definition depends linearly on the definition of the 
Heisenberg time. Now the final step: By inspection of 
the dynamics of the Floquet quantal system [Eqs. dl]),®] 
one can see (see also the derivation in the Stockmann's 
book [l| ) that the value of the variance of the angular mo- 
mentum at the point of stopping the diffusion ((Ap)^) is 
proportional to h^l'^, and to achieve equality we intro- 
duce a dimensionless numerical (empirical) factor by 
writing ((Ap)^) = h^l'^/a^, which on the other hand 
must be equal just to the classical value at stopping time FIG. 
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1. [Colour online] The phase portrait of the standard 

namely equal to {T/I)^Df,li^. From this it follows map [Eq. ©] for = 5 in (a) and = 7 in (b). In 

(b) we can see two tiny islands of stability of a period one 
orbit near {6, P) = (« 4.25, 0) and {6, P) = (« 4.25, 2-k) (see 
D (K^ \ tey± for more discussion). The other stable fixed point at 

hoc = ( ctj, ^ ^ I . (17) (61, P) = 2.13, 0) is hardly visible. 



immediately 



The numerical constant is found empirically by nu- 
merical calculations, for instance in the literature the 
case K = b with = 1 is found to be ai = 0.5 (how- 
ever, we find numerically ai = 0.45, taking into account 
Eq. ()17|) when studying the model's localization proper- 
ties). Thus, we have the theoretical formula for the local- 
ization length in the case of generalized classical diffusion 
[Eqs. P^ . (|13p ]. which we use in defining the scaling pa- 
rameter [Eq. (I5T]) ]. 

As for the classical system [Eq. ([3])], we mention that 
the fraction of the regular part of the classical phase space 
has been systematically explored using the GALI method 
[53] for the distinction between chaotic and regular clas- 
sical motion and its quantification for simple (and even 
for coupled) standard map(s) (see [53-l56| and references 
therein) , showing that this fraction decreases with K rel- 
atively slowly, then faster around Kcrit, has some smaller 
oscillations, at = 5 still amounts about 2.2%, and fi- 
nally for if > 7 it is zero for all practical purposes (much 
less than 1%). Thus, at if > 7 we have no problems with 
the effects of the divided classical phase space. However, 
there are important subtleties about the classical diffu- 
sion process and Df^{K) which we now discuss. 

We show the phase portraits of the standard map from 
Eq. ^ ioT K = 5 and for if = 7 in Fig. [TJ in order to 



demonstrate that at if = 5 we still have islands of sta- 
bility of relative area about 2.2%, which means that the 
effects of divided phase space cannot be neglected, whilst 
at if = 7 there are no large islands of stability. How- 
ever, there are still two tiny islands of stability around 
(6», P) = (« 4.25, 0) and (6*, P) = (w 4.25, 27r) (near a pe- 
riod one stable orbit) whose relative fraction in the phase 
space is found with the GALI method to be ~ 0.0162%. 
The if = 7 is the most widely used parameter value for 
the model in this paper. 

In case of the normal diffusion /i = 1 the theoretical 
value of Di{K) is given in the literature, e.g. in [sl] or 



Di{K) 



iif2 [1 - 2 J2(if) (1 - J2(if ))] , if if > 4.5 
0.15(if -ifcr)^ if is'cr < if < 4.5 



(18) 



where Kent — 0.9716 and J2{K) is the Bessel function. 
Here we neglect higher terms of order if~^. However, 
there are many important subtle details in the classical 
diffusion further discussed below. 

The dependence of the diffusion constant for the 
growth of the variance of the momentum on if is 
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very sensitive, and described in the theoretical result 
[Eq. ([T5)) ]. and fails around the accelerator mode intervals 
(27rn) < K < y/{2Tmy + 16, n any positive integer. In 
these intervals for the accelerator modes n = 1 we have 
two stable fixed points located at p = 0, 6 — n — Bq and 
p — 0, 9 = n + 6q, where 6*0 — arcsin(27r/if ). There 
are two unstable fixed points at p = 0, 6 = 9q and 
p = 0, 9 = 2Tr — 9o. In our case K = 7 oi Fig. [Hd we have 
6*0 — 1.114. Moreover, as the diffusion might even be 
anomalous, we have recalculated the diffusion constant 
numerically, and the results are shown in Fig. [2l We see 
that the dotted theoretical curve stemming from Eq. (fT8|) 
describes the diffusion constant well outside the acceler- 
ator mode intervals. In general, however, the diffusion 
might be non- normal, described in Eq. (jl2p . For the case 
K ~ 7, which is the main case that we study classi- 
cally and quantally in this paper, we find three different 
regimes of diffusion as shown in Fig. [31 As it will be seen 
below, for our purposes the middle regime with fi « 0.9 
and 169.82 is relevant and important, because loo 

is in the range of N, the dimension of the A^-dimensional 
matrix model introduced in the next section. 




IV. THE FLOQUET MODEL SYSTEM: A 
FINITE UNITARY MATRIX AS THE FLOQUET 
OPERATOR (IZRAILEV MODEL) 

The motion of the QKR [Eq. Q] after one period T of 
the ijj wavefunction can be described also by the following 
symmetrized Floquet mapping, describing the evolution 
of the kicked rotator from the middle of a free rotation 
over a kick to the middle of the next free rotation, as 
follows 

ij{9,t + T) ^Uij{e,t), (19) 

Thus, the 'il>{9,t) function is determined in the middle of 
the rotation, between two successive kicks. The evolu- 
tion operator U of the system corresponds to one period. 
Due to the instant action of the kick, this evolution can be 
written as the product of three non-commuting unitary 
operators, the first and third of which correspond to the 

free rotation during half a period G{t/2) — exp {ij^ffi^ , 

T = hT/I, while the second B{k) — exp{—ik cos 9) , 
k = Vq/H, describes the kick. The system's behavior 
depends only on two dimensionless parameters, namely 
T and fc, and its correspondence with the classical sys- 
tem is described by the relation K = kr = VqT/I. In 
the case K = kr ^ 1 the motion is well known to be 
strongly chaotic, for K > 7 almost without any regular 
islands of stability, as explained in the previous section. 
The transition to classical mechanics is described by the 
limit fc oo, T while K = const. We shall consider 
mostly the semiclassical regime k > K , where r < 1. 



° 5 10 15 20 25 30 35 

K 

FIG. 2. Diffusion in the Chirikov map. We show the value of 
the classical diffusion constant as a function of K, for two 
discrete times n, i.e. the number of the iterations of the 
standard map, n — 1000 in (a) and n = 5000 in (b). The 
smooth background (dotted) agrees perfectly with the theory 
[Eq. (|18p ]. whilst the peaks are due to the anomalous diffu- 
sion associated with the accelerator modes and other sticky 
objects around them. We have used 5000 initial conditions 
uniformly distributed in a grid of a square unit [0, 1] x [0, 1]. 
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Logl0(") 

FIG. 3. The variance of the momentum P in the standard 
map [Eq. ([Sjl] with K — 7 for the same initial conditions as in 
Fig. [2] as a function of the discrete time n (number of itera- 
tions), in log-log representation. The three slopes associated 
with different types of diffusion are /i = 1.7 (dotted), = 0.9 
(solid) and fj, — 1.45 (dashed). 
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In order to study how the locahzation affects the sta- 
tistical properties of the quasienergy spectra, we use the 
model's representation in the momentum space with a 
finite number TV of levels [H M, 



N 



U.n{t + T)^ J2UnrnU,n{t), n, m = 1, 2, iV . (20) 
m— 1 

The finite symmetric unitary matrix [/„,„ determines the 
evolution of an A^-dimensional vector, namely the Fourier 
transform Un{t) of ip{(^, t), and is composed in the follow- 
ing way 



(21) 



where Gw — exp is a diagonal matrix corre- 

sponding to free rotation during a half period T/2, and 
the matrix Bn'm' describing the one kick has the follow- 
ing form 



Br, 



1 



2iV-|-l' 



2N+1 

E 

1=1 



X exp 



{n! — m') 



27r/ 



2iV + 1 



(n' + m!) 



27r/ 



ik cos 



/ 2ttI 



V2A^+ 1 



2A^+ 1 
(22) 



The model in Eqs. ()20ll22p . which we refer to as Izrailev 
model, with a finite number of states is considered as the 
quantum analogue of the classical standard mapping on 
the torus with closed momentum p and phase 0, where 
Unm describes only the odd states of the systems, i.e. 
^{9) = —ijj{—6), provided we have the case of the quan- 
tum resonance, namely r = Aur / {2N + 1), where r is a 
positive integer, as in Eq. (|10l) . The matrix [Eq. (|22l) ] is 
obtained by starting the derivation from the odd-parity 
basis of sin(n0) rather than the general angular momen- 
tum basis exp(m6'). 

Nevertheless, we shall use this model for any value of 
r and k, as a model which in the resonant and in the 
generic case (irrational T/(47r)) corresponds to the clas- 
sical kicked rotator, and in the limit N ^ oo approaches 
the infinite dimensional model [Eq. (IT51) ]. restricted to 
the symmetry class of the odd eigenfunctions. It is of 
course just one of the possible discrete approximations 
to the continuous infinite dimensional model. 

The difference of behaviour between the generic case 
and the quantum resonance shows up only at very large 
times, which grow fast with (2A'^ + l)j as explained in 
section |lTl It turns out that also the eigenfunctions and 
the spectra of the eigenphases at finite dimension N of 
the matrices that we consider do not show any significant 
differences in structural behaviour for the rational or irra- 
tional T/(47r), which we have carefully checked. Indeed, 
although the eigenfunctions and the spectrum of the 
eigenphases exhibit sensitive dependence on the parame- 
ters T and k, their statistical properties are stable against 



the small changes of r and k. This is an advantage, as in- 
stead of using very large single matrices for the statistical 
analysis, we can take a large ensemble of smaller matrices 
for values of r and k around some central value of r = rg 
and k = ko, which greatly facilitates the numerical cal- 
culations and improves the statistical significance of our 
empirical results. Therefore our approach is physically 
mea ning ful. Similar approach was undertaken by Izrailev 
(see (Sg and references therein) . In Fig. |4] we show the 
examples of strongly exponentially localized eigenstates 
by plotting the natural logarithm of the probabilities 
Wn — \un\'^ versus the momentum quantum number n, for 
two different matrix dimensions N. By calculating the 
localization length loo from the slopes a of these eigen- 
functions using Eq. (|lip we can get the first quantitative 
empirical localization measure to be discussed and used 
later on. Here = 2/cr « 2.5 < iV = 398 for Fig. gji 
and « 2.2 < TV = 796 for Fig. H). 

At larger k the localization length can become compa- 
rable to N and the size effects start to play a role, there- 
fore loo is more difficult to determine as the fluctuations 
of 2/(7 are larger, but still can be done to some extent. 
Such a case is shown in Fig. [5^ with Zoo = 2/(t f« 57 <C 
N = 398 and in Fig. [5}) « 20 < = 796, where the 
quantification is difficult even if the localization length 
2/ct is well below N. We must be aware of the fact that 
the fluctuations in 2/ a are very large, as was observed al- 
ready in the pioneering works of Chirikov, Casati, Izrailev 
and Shepelyansky. 

Nevertheless, as long as the localization length is small 
enough, it is correctly predicted by the theory [Eq. ([T7|) ]. 
and is independent of the dimension N. As A^ increases 
with loo being flxed we approach the regime of strong 
localization and Poissonian statistics for the eigenvalues. 
Only when loo becomes comparable to A^ or larger, we 
observe gradual transition to the full quantum chaoticity, 
namely to GOE (or COE) behaviour, and anything in 
between. This transition is of central interest and is the 
main subject of the next sections. 

In Fig. we show the relationship between the 2/cr 
and the theoretical loo for variety of matrices [Eqs. ([^01 
[22]) ] and K = 7. We have statistically significant linear 
relationship with the slope w 0.9 close to unity (identity). 
Thus the theoretical loo from Eq. ([T7| agrees reasonably 
with the empirical localization length 2/cr as long as they 
are both sufficiently smaller than A^. 

We can define also the relative exponential localization 
measure, defined as the ratio of 2/ a and A^, 



A^ct' 



(23) 



which has meaning only if & <C 1, because 2/ a is well 
defined only if it is much smaller than A^, but even then 
we must be aware of large statistical fluctuations of 2/cr. 

There is another way to empirically quantify the de- 
gree of localization based on the information entropy, less 
sensitive to the finite size effects, denoted by j3\oc, and is 
discussed in the next section. In Fig. ^jp we show the 
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FIG. 4. (a) A sample of strong localized eigenstates for K — 7, 
r = 222, k ^ 2.00 and iV = 398 (b) Same ior K = 7,r = 444, 
k ^ 2.00 and N = 796. 



FIG. 5. (a) A sample of weak localized eigenstates for K = 7, 
r = 63, fc 7.05 and N = 398 (b) Same for /S' = 7, r = 127, 
k « 6.99 and iV = 796. 



relationship between b and Pioc- It is seen to be a linear 
relationship, but it must be emphasized that it applies 
only to sufficiently small values of b, because empirical 
exponential localization length 2/ a completely loses its 
meaning at values > A'^, or in fact, even much earlier. 
This is the reason why the slope in Fig. |61d <j ~ 0.404 
does not agree with Eq. (|45|) of section IVIH namely the 
slope there is I/7 « 0.25, 7 = 4.04, which is due to large 
fluctuations in a. 



V. DYNAMICAL (CHIRIKOV) LOCALIZATION 
OF THE EIGENSTATES AND ITS MEASURE 

Based on the examples of the eigenstates shown in 
the previous section, and following (36| and the refer- 
ences therein, we introduce the information entropy of 
the eigenstates as follows. 

For each A'^-dimensional eigenvector of the matrix Unm 
the information entropy is 



un) 



N 



where w„ 



^and E J 



Wn In Wn, 



1. 



(24) 



In case of the random matrix theory being applica- 
ble to our system [Eqs. dH]) and (|20II22I) ]. namely the 



COE (or GOE), due to the isotropic distribution of the 
eigenvectors of a COE of random matrices, we have the 
probability density function of on the interval [0, 1], 

■°""-"= .j«rw2) "-i""i'''"""'- <^^' 

It is easy to show that in the limit iV — > 00 this becomes 
a Gaussian distribution 



WN(\Un\) = Y ~^^P 



N\i 



(26) 



and the corresponding information entropy [Eq. (|24|) ] is 
equal to 



GOE 
N 



In ( iiVaj + 0{1/N), 
(27) 



where a — 



0.96, while ip is the digamma 



cxp(2-7) 

function and 7 the Euler constant (~ 0.57721...). For 
a uniform distribution over M states w„ = 1/M we get 
Jfjy ~ logM, and thus M « exp{J^]\[). Thus, we get the 
insight that the correct measure of localization must be 
proportional to exp(Jif]^), but properly normalized, such 
that in case of extendedness (GOE/COE) it is equal to 
N. 
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FIG. 6. 2/a vs loo, where a is the slope of the eigenfunctions 
in momentum space, namely of w„ = \unf vs n, for several 
it- values and iV = 398(n) and iV = 796(+) (fit slope « 0.89). 
(b) b vs. /3ioc, = 398(n) and iV = 796(+) (fit slope ^ 
0.404). 

Therefore the entropy localization length Ih is defined 

as 



Ih = A^exp - Jif^ 



GOE\ 
N ) 



(28) 



Indeed, for entirely extended eigenstates Ih — N . Thus, 
Ih can be calculated for every eigenstate individually. 
Hovifever, all eigenstates, while being quite different in 
detail, are exponentially localized, and thus statistically 
very similar. Therefore, in order to minimize the fluctu- 
ations one uses the mean localization length d = (Ih), 
which is computed by averaging the entropy over all 
eigenvectors of the same matrix (or even over an ensem- 
ble of similar matrices) 



d = (Ih) = iVexp {(J^n) - ^^°^) ■ 
The localization parameter /3ioc is then defined as 
d _ (Ih) 



Aoc 



N 



N 



(29) 



(30) 



Its relationship to b at small values /3ioc is shown in 
Fig. |6}d and discussed in the previous section. 

The parameter that determines the transition from 
weak to strong quantum chaos is neither the strength 
parameter k nor the localization length Zoo , but the ratio 



of the localization length Zoo to the size N of the system 
in momentum p 



Zoo ^ ^ 

N N 



(31) 



where Zoo ~ ^loc, the theoretical localization length, was 
derived in Eq. P7|) . A is the scaling parameter of the 
system. This is one of the main results of the present 
work as it incorporates normal diffusion /i = 1 and the 
anomalous diffusion /i 7^ 1. The relationship of A to /3ioc 
is discussed in section IVlIl 



VI. THE QUASIENERGY SPECTRUM AND ITS 
STATISTICAL PROPERTIES 



In this section we study the statistical properties of 
the spectrum of the eigenvalues corresponding to the 
(i) 

eigenstates m„ , labeled by j, of the Floquet operator 



Fnm introduced in Eq. 



namely 



,3 

'in ■> 



(32) 



where due to the unitarity of Fnm all eigenvalues \j must 
be on the complex unit circle, A^- = exjp{i<j)j). Thus, our 
analysis is concerned about the statistical properties of 
the spectrum of the eigenphases € [0, 2-k). This is 
equivalent to the quasienergies hwj = h(j>j/T |46| . where 
T is the period of the time-periodic Hamiltonian. As 
mentioned in section |TT1 the spectrum can be contin- 
uous (in case of the quantum resonance with rational 
T / (An) = r/q, having extended states, and no dynami- 
cal localization), or discrete (in the generic case of suf- 
ficiently irrational t/(Att) ) with dynamical localization. 
Even then, the spectrum has infinite level (eigenvalue) 
density as we have infinitely many eigenphases 0j on the 
interval [0,27r), and consequently all level spacings are 
zero. However, for any finite dimension N, no matter 
how large, the mean level spacing is 27t/N, and thus 
finite, and we can begin the statistical analysis of the 
quasienergy spectrum. 

Therefore, instead of using the infinite dimensional sys- 
tem [Eq. ([HI)], we study the finite dimensional Izrailev 
model [Eqs. (|20j|22p ]. It is a discrete approximation to 
the exact initial infinite dimensional system [Eq. ([S])], or 
in its symmetrized version [Eq. (I19p ]. In this case we can 
best use irrational t/{Att), to ensure that we are away 
from the quantum resonance. 

The most important statistical measure of the eigen- 
values is the level spacing distribution. We study only the 
level spacing distribution, but in three different ways of 
analyzing it, namely as the probability density function 
P{S), its cumulative distribution function W{S) and the 
so-called {/-function, introduced by Prosen and Robnik 
in [H,!!!] (see Appendix 
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In so doing we come to the central point of this work, 
namely the empirical almost identity relationship be- 
tween localization measure P\oc of the eigenstates and the 
spectral level repulsion parameter (exponent) in the level 
spacing distribution P{S), proposed already by Izrailev 
in [sll, [m and [11]. The localization length Zoo, derived 
by the semiclassical argument [Eq. ^T7\ ]. through the /3ioc 
- A relationship directly gives a prediction for the level 
repulsion parameter in the semiclassical regime oi k > K 
and large N. 

Here we reproduce all Izrailev's findings (for K = 5), 
generalize them (for many other values of K, predomi- 
nantly ai K = 7), sharpen his results and put them in 
the broader perspective including the autonomous (time- 
independent) Hamiltonian systems. We find that the 
Brody distribution captures the numerically found level 
spacing distribution statistically highly significantly, no- 
tably better than the distribution function proposed in 
^5, 36, 4i|, also [s^l, and this is entirely in line with the 
results on the dynamically localized chaotic eigenstates 
in time-independent Hamiltonian systems, like billiards 
psl m, nil , even in the mixed type regime (after sepa- 
rating the regular and chaotic eigenstates). 



assess the goodness of the fit of the theoretical level spac- 
ing distribution. 

The important special level spacing distributions that 
we are using are the following. 

• The Poisson distribution 

Ppoisso„(5) = e-^, l^Poisson(5) = l-e-^. (35) 

• The COE or GOE distribution 

PcOe{S) ~ Pwignor(S') = ^ CXp (^-^S"^^ , 

WcoEiS) = l-exp(^^^^. (36) 

• The Brody distribution [H, Hi] 



(37) 



where the two parameters Ci and C2 are deter- 
mined by the two generic normalization conditions 
that must be obeyed by any P{S), 



poo POO 

Level spacing distribution: P{S), W{S) and U{S) I p[S)dS = 1, {S) = / SP{S)dS = 1, (38) 

Ja Jo 



To study the eigenvalue statistics of quantum Flo- 
quet systems and quantum maps one considers the eigen- 
phases 0„ e [0, 27r) defined by A„ = e*"^". In such 
case the spectral unfolding procedure is very easy, as 
the mean level density is N/{2tt), i.e. the mean level 
spacing is 2tt/N. The histogram of the level spac- 
ing distribution P{S) is the distribution of the spacings 
Sn ■■= ^i(l>n+i - 4>n) with n = 1, N and (I>n+i ■= 0i, 
in the bins of certain suitable size AS'. The factor N/2tt 
ensures that the average of all spacings S^ is 1, and thus 
P{S) is supported on the interval [0,A^], and its upper 
limit goes to oo when — > oo. 

The cumulative distribution VF(5), or integrated level 
spacing distribution, preserving the full accuracy of all 
numerical eigenvalues/spacings, useful especially when 
the number of levels N is small, is defined as 



thus with (5*) = 1 being the mean distance between 
neighboring levels (after unfolding). Hence 



Ci = (/? + 1)C2, C2 



P + 2 



13+1 



(39) 



where T{x) denotes the Gamma function. In the 
strongly localized regime /3 = we observe Poisso- 
nian statistics while in the fully chaotic one f3 — 1 
and the RMT applies. The Brody cumulative level 
spacing distribution is 



W^br(S) = 1-cxp(-C2S'3+1). 



(40) 



W{S) = 



fW^. ,33, 



Finally, we shall use also the so-called [/-function (see 
the Appendix R)) 



U{S) — — arccos 

TT 



(34) 



• Izrailev distribution: In [35|, |4l|, Izrailev suggested 
the following distribution in order to describe the 
intermediate statistics, i.e. the non-integer /3 in the 
following PDF could be associated with the statis- 
tics of the quasienergy states with chaotic localized 
eigenfunctions, also approximating the level spac- 
ing distribution arising from the Dyson COE joint 
probability distribution [58] 



The {/-function has the advantage that its expected sta- 
tistical error SU is independent of S, being constant for 
each S and equal to SU = 1/ {-k^/N s)^ where A^^ is the 
total number of objects in the W{S) distribution. The 
numerical pre- factor 2/ti in Eq. (j34l) is determined in 
such a way that U{S) € [0, 1] when 1^(5') G [0, 1]. The 
{/-function is an excellent and refined criterion used to 



Piz{S) 



exp 



B - ^TT/? ) 5 



(41) 



where again the two parameters A and B are de- 
termined by the two normalization conditions (1) = 
(S) = 1 given above. 
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Of course, we must be fully aware of the fact that 
both, Brody and Izrailev distributions, are approxima- 
tions. It is clear that at /3br = 1 we get precisely 
Wigner surmise [Eq. which is the exact GOE only 

for two-dimensional Gaussian random matrices, and thus 
only an (excellent) approximation for the infinite dimen- 
sional GOE case. Indeed, if we try to fit the exact infi- 
nite dimensional GOE level spacing distribution with the 
Brody distribution, we do not get /3br = 1, but instead 
/?BR = 0.953, see [13] ■ Also, we should mention that 
Izrailev et al have published an improved distribution 
function [s^], which we have also tested, and is defined 
by 

Pi"r (5) = ASP{1 + Bpsy^^^ exp[-^/?52 - |(1 - ^S)] 

(42) 



where /(/3) 



/3 



0.16874 and A,B are the nor- 



malization parameters. We found (see below) that in our 
applications it is even worse than the original version 
(HH). 



Analysis of the level spacing distribution of 
numerical spectra 



B. 



For the numerical calculations and results regarding 
the spacing distributions P{S) (and W{S)) for the eigen- 
phases (t>j, we have considered a range of 41-values of the 
quantum parameter k (= 2, 3,..., 42) keeping fixed the 
classical parameter K — 7 (where the phase space is fully 
chaotic as shown in section|TTl. In order to ameliorate the 
statistics we considered a sample of 161 matrices Unm of 
size N = 398 (« 64,000 elements), in a similar manner 
as in [3^) with slightly different values of k (with the 
step size Ak — ±0.00125 <C k). For some samples we 
reached up to 641 matrices Unm of size N = 398 acquir- 
ing qualitatively the same results. 

For the ensemble of M = 641 matrices of size N = 398, 
in case K = 7 and k — 11, using the best fitting 
procedure (described in more detail below) we found 
/3br = 0.421 and all three representations clearly show 
excellent agreement with the best fitting Brody distri- 
bution. In Fig. [7^ we show the histogram. It is seen 
that Brody distribution is better fitting the data than 
the Izrailev distribution. Since the deviations are really 
small, the statistical significance very high, we plot in 
Fig-Eb the differences C/(data) — J7 (Brody /Izrailev) ver- 
sus Ubr- Thus if data are on the abscissa the agreement 
is perfect. As can be seen, the deviations from that are 
really small, and clearly smaller for Brody. In Fig. [T]: we 
show the fine differences of (data) — ly (Brody /Izrailev) 
versus W^br, and again we clearly see that Brody is bet- 
ter. Finally, in Fig. [7ji we show the same thing as in 
Fig. [Tj:, except not versus W but versus S instead. It 
must be emphasized that the improved Izrailev distribu- 
tion (I42p exhibits even larger deviations from the data 



than the original one (|4T|) . as depicted in Fig. [7] by the 
outer most line in panels (b,c,d). This is the main reason 
why we have not considered the new Izrailev distribution 
any further, and also in order to be compatible with 
the previous related results in the literature. 

More data for various k are shown in Table [D The in- 
dex PDF or CDF means that the fitting was done by 
using P{S) or W{S), respectively. The tick marks indi- 
cate which fitting is statistically better, based on the 
procedure. In order to illustrate some more cases from 
the Table HI namely for other values of /3, we show the 
results in Appendix [Bj Even more data are in the Table 

m 

Finally, we do a similar analysis for large matrices 
N = 4000, and take AI = 9 oi them in an ensemble. The 
results in Appendix |B] clearly show that Brody distribu- 
tion is an excellent fit to the level spacing distribution in 
all three representations, P{S), W{S) and U{S). 



C. Residues and test 

In the best fitting procedure we have calculated both 
the residues and the as follows: 

• PDFs residues 



N 



^PDFs = ^(pBR,iz(i) - data(i))" 



2 _ (-Pbr,iz(») - data(z))" 

XpDFs — / , 5 T^N 



(43) 



(44) 



In Fig. |5]we show three examples for the as a func- 
tion of the fitted parameter /3br or /?iz for the data from 
the Tables H UlIII It is clearly demonstrated that Brody fit 
is significantly better than Izrailev. 



VII. THE SCALING BEHAVIOUR OF Aoc 
VERSUS A 



It is well established that the degree of localization 
and the value of the spectral level repulsion parame- 
ter /3 are related. The parameter that determines the 
transition from weak to strong quantum chaos is neither 
the strength parameter k nor the localization length Zoo, 
but the ratio of the localization length Zoo to the size N 
of the system in momentum p, the scaling parameter A 
[Eq. (|3T|) ]. Here we present the following results 

• In plotting /3br versus /3ioc as shown in Fig. [9] for 
a great number of matrices at various parameter 
values we clearly see the linear relationship very 
close to identity. 



12 








Brody 






\ / 

w 


V>^X Izrailev 




D 2 4 6 


0.8 




(d) 


A 

Brody 




)L\r 








Izrailev 




D 1 2 3 4 


5 6 



FIG. 7. Intermediate statistics (panel (a)) for distribution P{S) (histogram - black line) of the model fitted with distribution 
Pbr (blue-solid line), Piz (red-dashed line pointed by the arrow in the inset figure) and P{^™ (black-dotted line) for M x N — 
641 X 398, K = 7 and = 11 (see text for discussion). In panels (b),(c),(d) we show the difference of the numerical data and 
the best fitting Brody (blue-thick line), "old" Izrailev (red- medium line) and Izrailev "new/improved" (black-thin line; always 
the outer one) PDFs by using the f/-function and Vt^-distribution. Thus, in case of the ideal fitting the data would lie on the 
abscissa. In this case, based on the P{S) fit we get /3br=0.421, /?iz(old)=0.416 and /3iz(new)=0.376, and based on the W{S) 
fit we get /3br=0.421, /3iz(old)=0.401, /3iz(new)=0.350. 
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FIG. 8. for Pbr and Pz for = 398 (a), iV = 796 (b) and = 4000 (c) for various values of the best fits for /3br (black 
fine) and /3iz (gray line) using 150 subintervals in 5*. 



• In Fig. [TU] we plot /3ioc versus A for the same en- 
semble of matrices as in Fig. [S] clearly, there is a 
functional relationship according to the scaling law 



Aoc — 



7A 



1 + 7A 



7 = 4.04, 



(45) 



which is similar to the scaling law [Eq. (^S)) ]. but 
not the same. Thus we see that when A — >■ cx) both 



/3ioc and /3br go to one, and we have extended eigen- 
states and GOE/COE spectral statistics, whilst in 
the limit A = we have strong localization, /3br 
and /3ioc are zero, and we have Poissonian spectral 
statistics. Fig. [S] shows what happens in between. 
The value 7 = 4.04 differs somewhat from 7 « 3.2 
in [5^, where /3ioc is plotted versus x « 4A. 

The quality of the fit of Fig. IH] is degraded a lot 
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TABLE I. Results of the best fitting procedure using Brody 
(BR) and Izrailev (IZ) distribution. The index PDF or CDF 
means that the fitting was done by using P{S) or W{S). The 
tick marks indicate which fitting (BR or IZ) is statistically 
better, based on the procedure. 



if = 7 and M X iV = 161 X 398 



k 


oPDF 

Pbr 


aCDF 

Pbr 


oPDF 

Piz 


oCDF 

Piz 


5 


0.131197 


0.139808 


0.121197/ 


0.109808/ 


8 


0.280887/ 


0.286922/ 


0.265887 


0.256920 


11 


0.421398/ 


0.420996/ 


0.416398 


0.400996 


14 


0.581301/ 


0.574325/ 


0.596301 


0.574325 


17 


0.670686/ 


0.659891/ 


0.705686 


0.679891 


20 


0.713984/ 


0.699906/ 


0.763984 


0.739906 


23 


0.791288/ 


0.778780/ 


0.861288 


0.838780 


26 


0.812993/ 


0.797852/ 


0.892993 


0.867852 


29 


0.832067/ 


0.821543 


0.912067 


0.891543/ 



TABLE IV. We show the values of fi, for the classical 
diffusion. The coefficients are needed for the quantum 
kicked rotator in estimating Zoo , and they are also determined 
numerically by seeking the best agreement in Fig. 1101 



K 


A* 






K 
<J 


n QQ 




13.182 


7 


« 0.90 


0.20 


169.82 


10 


« 1.00 


0.50 


31.68 


12 


« 1.00 


0.50 


87.09 


14 


« 1.00 


0.50 


134.89 


17 


« 1.01 


0.50 


97.05 


20 


« 0.99 


0.50 


275.42 


25 


« 0.98 


0.50 


405.50 


30 


« 1.00 


0.50 


389.04 


35 


^ 1.00 


0.50 


467.73 



TABLE II. Results of the best fitting procedure. The same 
as in Table [H but for A'^ = 796 and other parameter values k. 





K = 


7 and M x N 


= 81 X 796 




k 


oPDF 

Pbr 


aCDF 

Pbr 


oPDF 

Piz 


oCDF 

Piz 


6 


0.094946 


0.102980 


0.084946/ 


0.082980/ 


7 


0.122990 


0.135550 


0.11299/ 


0.1079550/ 


12.5 


0.316602/ 


0.319169/ 


0.306602 


0.289169 


15 


0.407302/ 


0.400864/ 


0.402302 


0.380864 


17.5 


0.500219/ 


0.491786/ 


0.505219 


0.481786 


20 


0.570621/ 


0.560884/ 


0.585621 


0.560884 


25 


0.696955/ 


0.680933/ 


0.741955 


0.710933 


30 


0.778601/ 


0.761853/ 


0.848601 


0.821853 


35 


0.831750/ 


0.816494/ 


0.911750 


0.886494 



when the size of the ensemble is decreased, as we 
have observed. 

• The quality of the fit when using /3iz instead of 
/3br is decreased, as we have checked carefully and 
as shown in Tables U [TT] and IIIIl 



TABLE III. Results of the best fitting procedure. The same 
as in Table [II but for A'^ = 4000 and other parameter values 
k. 



K ^7 and M xN = 9 x 4000 



k 


oPDF 

Pbr 


aCDF 

Pbr 


oPDF 

Piz 


oCDF 

Piz 


10 


0.058254 


0.064513 


0.053254/ 


0.055413/ 


15 


0.123504 


0.136268 


0.113504/ 


0.116268/ 


20 


0.223351/ 


0.228405/ 


0.208351 


0.198405 


25 


0.303080/ 


0.304349/ 


0.293080 


0.274349 


30 


0.424033/ 


0.410826/ 


0.419033 


0.390826 


35 


0.497268/ 


0.485010/ 


0.502268 


0.475010 


40 


0.569422/ 


0.553461/ 


0.586422 


0.553461 


45 


0.614519/ 


0.596811/ 


0.639519 


0.606811 


50 


0.668207/ 


0.651155/ 


0.708207 


0.671155 


55 


0.694437/ 


0.679703/ 


0.739437 


0.709703 


60 


0.761739/ 


0.740879/ 


0.821739 


0.790879 



The numerical factor is determined by numerical 
calculations, namely by seeking the best agreement in 
Fig.[Tni It is interesting to note that in the case K — 5 ex- 
tensively studied by Izrailev ([36'] and references therein) 
we find /i = 0.99, which is compatible with /i = 1, 
but = 0.45 which is only approximately his value 
0.5. Then Di sa and K^t'^ = fc^ and we have 

Zoo ~ fc^/4. The data for all cases that we treat in this 
paper are given in Table ITVl 

In [5§| the following scaling law was proposed 

Pxoc — , where x= — and 7 « 0.8. (46) 

1+70; N 

A banded random matrix model has been proposed 
[37I l60l - [6^ to explain the above scaling relationship, 
based on the fact that Eq. (|5]) is an infinite banded matrix 
of band width approximately fc, which can be reduced 
to the finite dimensional model with the same property. 
However, one should observe the fact that k"^ /{AN) is 
just an approximate value of A valid for the special case 
when ^ — 1 (normal diffusion), ai = 1/2 and K = 
where Di « K'^/{2t'^) — In fact, we find numeri- 

cally ai — 0.45, not 0.5. Moreover, we must be aware of 
the fact that at ii' = 5 we still have islands of stability 
in the classical phase space, implying problems with the 
divided phase space in the quantum picture. This is the 
reason why we consider the cases K > 7, but nevertheless 
check Izrailev's results [s^ limited to if = 5 (Table ITV]). 
As seen in Fig. [31 as an example, in case K — 7 we have 
three different diffusion regimes, and for our purposes the 
middle regime with fi w 0.9 and « 169.82 is relevant 
and important, because Zoo is in the range of N. In the 
present work we thus generalize the empirical scaling law 
[Eq. ([45])] to all values oi K > Kent ~ 0.97 and show 
that the scaling /3ioc versus A persists under the proper 
generalizations, now including anomalous diffusion. 

Thus the knowledge of the theoretical scaling parame- 
ter A from Eq. (|3ip . in the semiclassical limit, enables us 
to calculate the spectral level repulsion parameter /3br 
(or /3iz) according to the scaling law [Eq. This law 

clearly shows that at fixed parameters K , t and k = K/t 
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FIG. 9. [Colour online] The fit parameter /3br as a function 
of Ploc for 161 X 398 elements for various values of K and for 
a wide range of k values. The best fitting straight line is very 
close to identity. 
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FIG. 10. [Colour online] The parameter /?ioc vs. A for 
various K and k, (the same as in Fig. O where the scaling 
law [Eq. (|45|) ] is shown with the black line. 



in the limit iV — ^ oo we always get /3br = 0, i.e. Pois- 
son distribution. On the other hand, when /qo becomes 
greater than N we see approach to the random matrix 
(GOE or COE) behaviour. The scaling properties of this 
section are one of the main conclusions of this paper. 
Recently, strong evidence has been found that a sim- 
ilar relationship can exist for the analogous quantities 
in dynamically localized time-independent Hamiltonian 
systems [28^291. 



VIII. DISCUSSION AND CONCLUSIONS 

In this paper we study in detail the relationship be- 
tween the localization measure of the eigenstates and 
the spectral level repulsion parameter in the quantum 
kicked rotator. First we confirm and improve the results 
of references [H, HI, [39l - l4l| for K = 5, and then we go 
substantially beyond his work by doing the analysis for 
many different values of if > 5 (dimensionless classical 



kick parameter), not only K = 5, namely K € [5,35], 
with various classical dynamical regimes, and many dif- 
ferent k € [5, 60] (dimensionless quantum kick parame- 
ter). Namely, we include also the cases with accelerator 
modes and generalized (anomalous) diffusion (subdiffu- 
sion and superdifFusion). 

The classical kicked rotator [Eq. ([IJ] is one of the most 
important model systems in classical and quantum chaos. 
We have studied the semiclassical regime where k > K 
and analyzed the eigenstates and eigenvalues (quasiener- 
gies, or eigenphases G [0,27r)), in particular the as- 
pects of the dynamical localization. The infinite dimen- 
sional case has a finite localization length loo (in the 
space of the angular momentum quantum numbers) , and 
exhibits the Poissonian statistics of the level spacings. 
We study a finite iV-dimensional model [Eqs. (|20ll22p ] 
proposed by Izrailev [s^ [H, Isol - lill ] for the case of odd- 
parity eigenstates. In this model the intermediate spec- 
tral statistics is observed, ranging from Poissonian to the 
RMT statistics (GOE/COE statistics), depending on the 
value of the scaling parameter A = loo/N. We have 
shown that A can be calculated theoretically in terms of 
generalized classical diffusion properties of the standard 

map expressed as Zoo = in Eq. (HZl), which 

is a major result of this paper as the anomalous diffusion 
{n 1) is now included. The degree of localization is also 
estimated by means of information entropy described by 
the parameter /3ioc, which goes from (Poisson) to 1 
(GOE), and is uniquely determined by A, as described 
by the scaling law /3ioc = t^^, 7 = 4.04 in Eq. p5)) . 

We find that in almost all cases the Brody distribution 
correctly captures the level spacing distribution at all 
values of the corresponding level repulsion exponent /3br, 
noticeably better than the distribution function proposed 
in [35_,,41j, and also in 37|, as demonstrated in section lVTl 
We confirm and significantly refine the result that /3br is 
identical to /3ioc within the statistical fiuctuations. 

These results have been obtained for the time- 
dependent system, the quantum kicked rotator, but we 
have evidence that similar conclusions can be reached 
in time-independent chaotic Hamiltonian systems, either 
in the mixed type regime or fully chaotic regime, when 
studying (after separation of regular and chaotic eigen- 
states) the localized chaotic eigenstates and the statistics 
of the corresponding chaotic (irregular) level sequences. 
This has been confirmed in the case of the billiard with 
mixed type dynamics 6^ [6^ in the recent paper by 
Batistic and Robnik [28 1 indirectly, and directly very re- 
cently [2^, [s^j. Moreover, it has been shown that differ- 
ent but equivalent localization measures can be intro- 
duced which are simply related to the Brody param- 
eter (Batistic and Robnik 2013). In case of billiards 
the quality of the spectral statistics is even much bet- 
ter than in the case of the quantum kicked rotator, due 
to the possibility to calculate a much larger number of 
high-lying eigenstates. Another important paradigmatic 
model is the hydrogen atom in strong magnetic field lei- 
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|68|, which in addition to the various bilhards is a good 
candidate for further theoretical and experimental stud- 
ies. 

To derive the fractional power law level repulsion and 
the emerging Brody distribution, as a consequence of the 
dynamical localization in chaotic eigenstates, is an im- 
portant open theoretical question. 
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we have 



Pik) 



k\{N, - k)l 



W''{l-W) 



(Al) 



Then the average values are equal to < k >= NgW, < 
k^ >^ NsW + N,{Ns-l)W^ and the variance ^(fc) =< 
P > - <k >2= N,W{1 - W). But the probability W 
is estimated in the mean as k/Ng. Its variance is 



V{W) = V 



and therefore the estimated error of W (standard devia- 
tion, the square root of the variance) is given by 



IW{1-W) 



Ns 

Transforming now from W{S) to 



(A3) 



(A4) 



/lv7 



U{S) = - arccos x/l - W{S), 

TT 

we show in a straightforward manner that SU 

and is indeed independent of S. From the (choice of the 
constant pre-factor in the) definition Eq. (jA4l) one sees 
that both U{S) and W{S) go from to 1 as 5 goes from 
to infinity. 



Appendix B: More spectral statistical analysis 



Appendix A: The U-function representation of the 
level spacing distribution 

First we estimate the expected fluctuation (error) of 
the cumulative (integrated) level spacing distribution 
W{S), which contains Ng objects. At a certain S we 
have the probability W that a level is in the interval 
[0, W] and 1 — W that it is in the interval [W, 1]. Assum- 
ing binomial probability distribution P{k) of having k 
levels in the first and Ng — k levels in the second interval 



In continuation of the section rVTl we show some more 
examples of the statistical spectral analysis for various 
values of /3br, in Figs. [TTl [T^] and [T21 Finally, we do a 
similar analysis for large matrices N — 4000, and take 
M — 9 oi them in an ensemble. The results are shown 
in Fig. [HI and more data are collected in Table Hill which 
are partially presented in . We clearly see that Brody 
distribution is an excellent fit to the level spacing dis- 
tribution in all three representations, P{S), W{S) and 
U{S). 
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